Chiral symmetry restoration and the Z3 sectors of QCD 
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Quenched SU(3) lattice gauge theory shows three phase transitions, namely the chiral, the decon- 
finement and the Z3 phase transition. Knowing whether or not the chiral and the deconfinement 
phase transition occur at the same temperature for all Z3 sectors could be crucial to understand the 
underlying microscopic dynamics. We use the existence of a gap in the Dirac spectrum as an order 
parameter for the restoration of chiral symmetry. We find that the spectral gap opens up at the 
same critical temperature in all Z3 sectors in contrast to earlier claims in the literature. 
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One of the most remarkable features of the QCD phase 
transition is the simultaneous vanishing of confinement 
and the restoration of chiral symmetry. Although there 
is much debate about the underlying mechanism which 
links the two transitions, their coincidence at a common 
critical temperature is considered a well established fact 
in lattice simulations (see e.g. ^ for a review). Interest- 
ingly, QCD in the quenched approximation has an addi- 
tional Z3 symmetry which allows us to obtain additional 
information. For pure gauge theory the phase transition 
can be described as spontaneous breaking of the Z3 sym- 
metry . For temperatures larger than Tc the Polyakov 
loop acquires a non-vanishing expectation value with a 
phase in one of the three sectors. Knowing whether the 
observed strict correlation between the chiral and decon- 
finement phase transitions persists in all Z3 sectors could 
provide important clues for the understanding of confine- 
ment. 

In it was claimed on the basis of lattice calculations 
with staggered fermions, that the restoration of chiral 
symmetry happens at different temperatures in the real 
{ip ^ 0) and complex sector of the Polyakov loop {(p ^ 
±2n/3). Several subsequent articles |Q| analyzed possible 
mechanisms to explain this observation. 

In this letter we reexamine this problem within lattice 
QCD using chirally improved fermions. In contrast to ||] 
we find that the critical P of the chiral phase transition 
does not depend on the Z3 sector, but is coincident with 
the Z3 breaking transition in all three sectors. 

Let us note that even if it were irrelevant for full 
QCD the center symmetry could still lead to fascinat- 
ing phenomenological consequences in supersymmetric 
Yang-Mills theories Thus the investigation of its 

properties is also relevant on more general grounds. 

Instead of directly measuring the chiral condensate we 
analyze in detail the spectrum of the lattice Dirac oper- 
ator. The density p(A) of the eigenvalues A of the Dirac 
operator is connected to the chiral condensate via the 
Banks-Casher formula [01, 



(1) 



where p(0) is the eigenvalue density near the origin and 



V is the volume of the box. Note that exact zero-modes 
which come from isolated instantons do not contribute 
to the density p(0) at the origin. The reason is that 
the number of zero-modes is believed to scale as 1/-'^/^ 
and thus they do not contribute when performing the 
thermodynamic limit in Eq. ([T]). At low temperatures 
when QCD is in the chirally broken phase the density 
is non-zero at the origin, while in the high temperature 
phase p is zero in a finite region around the origin, i.e. the 
spectrum develops a gap (up to isolated zero modes) and 
the chiral condensate vanishes (compare Fig. |^ below). 
The question whether chiral symmetry is restored at the 
same critical temperature in all sectors of the Polyakov 
loop can now be reformulated in terms of the spectral 
gap: As we increase the temperature, does the gap open 
up at the same temperature for all three sectors of the 
Polyakov loop? 

Before the advent of chirally symmetric formulations 
for the lattice Dirac operator such a study was quite awk- 
ward. In particular the spectrum of the Wilson lattice 
Dirac operator shows large fluctuations close to the ori- 
gin and the notion of a spectral density is not well de- 
fined. The situation has changed, since the re-discovery 
of the Ginsparg- Wilson equation |^ . Dirac operators D 
which obey the Ginsparg- Wilson equation have eigenval- 
ues which lie on a circle and it is straightforward to iden- 
tify a spectral density and study the emergence of the 
spectral gap. However, the only exact solution of the 
Ginsparg- Wilson equation, the overlap operator ||], has 
the drawback of being very expensive in a numerical im- 
plementation. 

Here we work with the chirally improved operator 
which is a systematic expansion of a solution of the 
Ginsparg- Wilson equation [ pO[ . In particular we use an 
approximation which has 19 terms in the expansion and 
is described in detail in ||ll[. The computation of the 
eigenvalues of the Dirac operator was done with the im- 
plicitly restarted Arnoldi method |p^ . 

For our quenched gauge configurations we use the 
Liischer-Weisz action ||l^. We work on lattices of size 
Lt X with the temporal extent Lt — 6 and two values 
for the spatial extent, L = 16 and L ~ 20. We use peri- 



1 



odic boundary conditions for the gauge fields, while for 
the fermions the boundary conditions are periodic only 
for the space directions but anti-periodic for the time di- 
rection. Our statistics is 800 configurations for the 6 x 16'^ 
lattices and 400 for the 6 x 20"^ lattices. We use 6 dif- 
ferent values of the inverse coupling /3 which gives rise 
to ensembles on both sides of the phase transition. In 
Table | we list our values of (3, the lattice spacing a 
and the temperature T. We used the coefficients given by 
tadpole-improved perturbation theory 15 1. Our val- 
ues for the couplings P^t and f3pg for the rectangle and 
parallelogram terms in the Liischer-Weisz action can be 
found in ITl. 



p 


8.10 


8.20 


8.25 


8.30 


8.45 


8.60 


a [fm] 


0.125 


0.115 


0.110 


0.106 


0.094 


0.084 


T [MeV] 


264 


287 


299 


311 


350 


391 



TABLE I. Parameters for our gauge field configurations. 
We list the values of /3, the lattice spacing a and the temper- 
ature T. 

Let us begin the discussion of the spectral gap with a 
look at typical spectra of our Dirac operator. In Fig. |l| 
we show the distribution in the complex plane of the 50 
smallest eigenvalues A for three different gauge configura- 
tions on 6 X 20"^ lattices. The symbols are our numerical 
results and the full curve is the so-called Ginsparg- Wilson 
circle, i.e. the circle of radius 1 in the complex plane with 
center 1. For our approximate Ginsparg- Wilson opera- 
tor the eigenvalues do not fall exactly on the circle but 
show small fluctuations around the circle. However, the 
eigenvalues are sufficiently well ordered to allow for the 
notion of a spectral density and a clear identification of 
the spectral gap. 

The plot on the left-hand side shows the spectrum for 
a configuration in the low temperature, chirally broken 
phase. For this case the eigenvalues extend all the way to 
the origin and there is a non- vanishing p(0) such that the 
Banks-Casher relation (||) gives rise to a non-vanishing 
chiral condensate. 

The central and the right-hand side plots show spec- 
tra for configurations in the high-temperature, chirally 
symmetric phase. The central plot is for a configuration 
with complex Polyakov loop P, while the right-hand side 
result is for real Polyakov loop. Both of these plots have 
a well pronounced spectral gap. The spectral density at 
the origin vanishes and so does the chiral condensate. For 
the complex sector the gap is considerably smaller than 
for the real sector. 

One can understand this difference between the sectors 
by considering the fcrmion boundary conditions. In the 
real sector the boundary condition 

V'(x,t+1/T) = -^{x,t) (2) 

gives a Matsubara frequency ttT to the fermions. In the 
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FIG. 1. Three typical spectra of the chirally improved op- 
erator. Only the 50 eigenvalues closest to the origin are plot- 
ted: Left plot: the chirally broken phase (6 x 20^, /3 = 8.10); 
Central plot: the symmetric phase (6 x 20^, f3 = 8.60) for com- 
plex Polyakov loop; Right plot: the symmetric phase for real 
Polyakov loop. The full curve is the Ginsparg- Wilson circle. 



complex Z3 sectors the boundary condition is effectively 
^ix,t + l/T) = e±*"/V(^,0, (3) 

giving a Matsubara frequency ttT/S. In the free-field 
case (i.e. (3 — > 00) the smallest eigenvalue is equal to 
the Matsubara frequency, giving a gap 3 times larger in 
the case of (H) compared with (g). It is thus reasonable 
that the real sector gap is considerably larger than the 
complex sector gap in the interacting case too. 

In quenched QCD the finite temperature phase transi- 
tion appears to be a weak first order phase transition . 
A first order phase transition is governed by the mixing 
of two phases and the behavior of their free energies. In 
our particular example we have a low temperature phase 
characterized by a vanishing spectral gap and a high tem- 
perature phase with a finite spectral gap. For tempera- 
tures sufficiently below or above the critical temperature 
the system is in only one of the two phases while near 
the critical temperature the system shows mixing of the 
two phases. 

We demonstrate this mixing in Fig. ^ where we show 
histograms for the distribution of the spectral gap at 
three different values of the temperature. We define the 
spectral gap g\ to be the imaginary part of the smallest 
eigenvalue which is not a zero-mode (as remarked above, 
zero-modes do not contribute to the spectral density). 
Zero-modes can be identified uniquely since for our chi- 
rally improved operator it can be shown that they have 
exactly vanishing imaginary part and the corresponding 
eigenstates have a non- vanishing matrix element with 75 , 
while this matrix element vanishes identically for non- 
zero-modes. 

We show histograms for the distribution of g\ for 
T < Tc,T - Tc and for T > T^. The top row displays 
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FIG. 2. Histograms for the spectral gap ag\. We show our 
results for lattice size 6 x 20^ at three temperatures T < Tc 
{(3 = 8.10), r ~ Tc {P = 8.30) and T > Tc (P = 8.60). 
The top row displays the results for the complex sector of the 
Polyakov loop P, while the bottom row is for real Polyakov 
loop. 



the results for the complex sector of the Polyakov loop 
P, while the bottom row is for real Polkyakov loop. The 
data were computed on 6 x 20'^ lattices. Since for the real 
sector {(fi ~ 0) the statistics is only half of the statistics 
for the complex sector (i^ ^ +2tt/3 and —2tt/3) we dou- 
bled the bin size for the two histograms in the real sector 
at T - Tc and T > T^. 

At T < Tc we find for all sectors a single peak near the 
origin. This peak is not located exactly at since also in 
the chirally broken phase the Dirac operator of a finite 
system has a microscopical gap which vanishes as 
p7| . For temperatures near Tj, the histograms show a 
clear double peak structure characteristic for the first or- 
der transition. The left peak corresponds to the chirally 
broken phase with a vanishing gap and the right peak 
is from the chirally symmetric phase with non-vanishing 
gap. As one increases the temperature further, only the 
right peak survives. As already noted in the discussion 
of Fig. |l| the gap is larger in the real sector, i.e. the right 
peak sits at larger values of g\ for the real sector. In ad- 
dition this peak is wider than the corresponding peak in 
the complex sector, i.e. the gap fluctuates more strongly 
around its mean value in the real sector. 

In order to describe the first order transition we use a 
simple ansatz for the behavior of observables. Let us first 
discuss the somewhat simpler case of the Polyakov loop. 
In an infinite system the Polyakov loop P vanishes below 
f3c and has a non- vanishing modulus |P| above f3c- On a 
finite lattice the Polyakov loop does not vanish exactly 
below Pc but disappears like V~^^^, i.e. like cL^^^^ with 
some constant c. Above /3c the dependence on L is negli- 
gible and to leading order P is linear in /3, i.e. described 
by d -I- k{(3 — (3c)- Following the ideas in [|l8[ one arrives 
at the conclusion that near the transition the expecta- 
tion value of the modulus of the Polyakov loop should be 
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FIG. 3. The expectation value (|P|) of the modulus of the 
Polyakov loop as a function of /3. The symbols indicate the 
numerical results while the full curve is a fit to Formula (^. 
We display results for 6 x 16^ and 6 x 20^ lattices. 



given by: 

_ cL-^/^e-^f^'<^0-P^) +i[d+k{P~Pc)] 

The term —AfV{f] — f]c) is the difference in the free ener- 
gies of the two phases. At (3 — (3c the two free energies are 
equal while below Pc the free energy of the chirally bro- 
ken phase is smaller than the free energy of the chirally 
symmetric phase and vice- versa above /3c ■ The factors 3 
in the second terms in the numerator and denominator 
come from the three possible values for the phase of the 
Polyakov loop. In Fig. ^ we show a fit of Formula (Q) 
(full curves) to our numerical data (symbols). In partic- 
ular we present a common fit to both the 6 x 16"^ and 
6 X 20'^ ensembles. This is possible since the parameters 
d and k are essentially independent of L. The fit result 
for /5c is given in Table |^ below. The fit demonstrates 
that both the L and the (3 dependence are well described 
by Formula (^). 

For the expectation value of the spectral gap 17^ we use 
a similar ansatz. 



c'L- 



(5) 



-A/L^(/3-/3e) 



i[drAL) + krAL){(3-(3c)] 



The subscripts r and c indicate the real respectively com- 
plex sectors of the Polyakov loop. Note that now we use 
the known L^^ behavior of the microscopical spectral gap 
in the chirally broken phase [ p^f . In the chirally symmet- 
ric phase the gap is essentially linear in (3 but the coeffi- 
cients dr^c and kr^c turn out to be L-dependent. Thus in 
a common fit to the 6 x IQ^ and 6 x 20^^ ensembles these 
parameters had to be varied independently. Again the 
fit results for (3c are given in Table ||. In Fig. W we plot 
our numerical data for the gap together with the curves 
(^) . The top plot gives the results for the complex sector 
while the bottom plot shows the real sector. As for the 
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FIG. 4. The spectral gap in lattice units as a function of 
p. The symbols indicate the numerical results and the full 
curve is a fit to Formula (P). We display results for 6 x 16'^ 
and 6 X 20^. 



Polyakov loop we find that the numerical data are rea- 
sonably well described by the simple first order transition 
formula. 

When comparing the results for the critical beta as 
given m Table || we find that within the accuracy we 
achieved the spectral gap vanishes at the same (3c for both 
the real and the complex sectors of the Polyakov loop. 
Furthermore this value is compatible with (3c as obtained 
from the analysis of the Polyakov loop. Combining the 
three methods wc find a critical temperature of 300 ± 3 
MeV for the Liischer-Weisz action which is slightly larger 
than the result for Wilson's gauge action. 



Measurement : 




) complex 


{gx)real 


(3. ■■ 


8.24(1) 


8.29(2) 


8.27(2) 


n [MeV] : 


296(3) 


308(5) 


303(5) 



TABLE II. Values of /3c and the critical temperature from 
the analysis of the Polyakov loop and the spectral gap in 
the complex and the real sectors for the ensembles with 
Liischer-Weisz action and chirally improved fermions. 
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